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Abstract 

The correlation functions for a strongly correlated exactly solvable one-dimensional bo- 
son system on a finite chain as well as in the thermodynamic limit are calculated explicitly. 
This system which we call the phase model is the strong coupling limit of the integrable 
q-boson hopping model. The results are presented as determinants. 

Introduction 

In this paper we investigate a one-dimensional strongly correlated boson lattice system where 
the kinetic term depends on the occupation of the adjacent lattice sites. It was introduced 
and solved by the quantum inverse scattering method (QISM) in [jy, g, |[ [|. The natural 
dynamical variables for this model are the so-called q-bosons || closely related to the quantum 
algebra formalism [p[ . This model is called the q-boson hopping model. It can be treated 
as a quantization jj of the classical Ablowitz-Ladik equation || which is one of the possible 
integrable lattice versions of the nonlinear Schrodinger equation In the scaling continuous 
limit the g-boson hopping model becomes the Bose gas model. 

We calculate correlation functions in the special case of the lattice g-boson model correspond- 
ing to the infinite value of the coupling constant. It appears that this special case is related 
to the phase operators of the quantum nonlinear optics so we will call it the phase model [|o| . 
Usually the infinite coupling limit corresponds to the free fermion point (the most well-known 
example is the Bose gas model). On the contrary the phase model cannot be regarded as a free 
fermion model. 

The description of the correlation functions for the models solved by means of Bethe Ansatz 
is based on the representation for the correlators as the Fredholm determinants of linear integral 
operators. Such representations were obtained for the first time in jll], [l2| for the simplest 
two point equal-time correlators of one-dimensional impenetrable bosons. Later they were gen- 
eralized for the case of time-dependent correlators for the models which are the free-fermion 
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points of models solved by means of the quantum inverse scattering method (the Bose gas in 
the infinite coupling limit [Ojand the isotropic XXO Heisenberg chain |l4|, |l5| ]) and also for 
the case of finite interaction (Jl6|, [lTj] for the one-dimensional Bose gas and ]1§|] for XXZ chain) . 
Such representations give an opportunity to get classical integrable equations for correlators 
which can be used, in particular, to calculate the long time and large distance asymptotics of 
correlation functions @ |o[ §lj ||]. 

The main result of this paper is the determinant representations for correlation functions. 
We consider the finite chain as well as the thermodynamic limit at zero and finite tempera- 



ture. The form factors for the phase model were calculated in |23 . Some of the results in the 
thermodynamic limit were reported in fLOf . 

This paper is organized as follows. In the first Section the g-boson model is defined and 
different limits are discussed. In Section 2 the Bethe Ansatz solution for the g-boson hopping 
model is presented. In Section 3 we discuss the solution of the phase model and its thermody- 
namics. In Section 4 the scalar products are calculated. In Section 5 we derive the determinant 
representation for the simplest correlation function which is the "darkness formation probabil- 
ity" . In Section 6 we calculate the two-point time-dependent correlation functions for the finite 
lattice while the thermodynamic limit is obtained in the last Section. The derivation of the form 
factors is presented in Appendix. 



1 Q-Boson hopping model 

The q-boson hopping model is defined by the Hamiltonian 

1 M 

H « = -2 E( B » B «+1 + B " B n+l - 2N n), (1-1) 
n=l 

with the periodic boundary conditions, M + 1 = 1. Operators B n , Z?t , and N n — form the 
q-boson algebra: 

[Ni, B)] = B\Sij , [N h Bj] = -BiSij, (1.2) 
[B i ,B}]=q- 2N *S ij . 

The c-number q is taken to be q — e 1 . We consider real 7 > 0. This algebra was investigated in 
H and the connectio n wi th the quantum SU q {2) algebra was established in ||. 

The Hamiltonian (l.l)of the q— boson hopping model commutes with the operator of total 
number of particles, 

M 

N = J2N n , [H q ,N] =0. (1.3) 

71=1 

The Hamiltonian with the chemical potential /2 

H h =H q -flN, (1.4) 
which has the same eigenstates will be considered further. 
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The g— boson algebra (1.2) possesses the representation in the local Fock spaces formed by 
the g-boson normalized states 



B | 0)j = 0, Nj | 0), = 0, 
B] I nj}j = [nj + 1]* | rij + 1), , B 6 \ nj ) 3 



(1.5) 



'J I '"J/J ~ l'"J r A J I J 1 A /J > J I "J/J 

and the normalized states on the whole lattice are 

M 



I n j - 



(1.6) 



4=1 
M 



in>=ni n ^=nw)- i (5t) r 

4=1 4=1 
The notation [n]! = rifc=iM ^ s use d and the 'box' is 



0). 



1-9" 



II = 



i=i 



11 1 -q- 2 ' 

The q— boson operators ([l]^) can be expressed in terms of the ordinary canonical bosons 

b},b j ,[b i ,b]] = 6 ij ,N j = b}b j , (1.7) 



as 



Bj = (BW 



N j + 1 



(1.8) 



If q — > 1 (7 — > Q), the q— bosons become ordinary bosons 



B, 



At 7 = (the free boson limit) the Hamiltonian (1.1) becomes the linear hopping model Hamil- 
tonian 



M 



(1.9) 



while the g— boson algebra is the ordinary boson algebra. 

In the continuum scaling limit the lattice constant S — > 0, MS = L, 7 = ^C(5, and the g-boson 
Hamiltonian (1.1) is just the Hamiltonian of the Bose gas model with the repulsive S- function 
interaction of strength c : 



H= / dx{d x b f (x)d x b(x) + cb\x)tf (x)b(x)b(x)} , 
[b(x),tf(y)]=5(x-y). 
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Representation (1.8) shows that (1.1) involves non-linearities of all orders about the ordinary 
boson hopping model (1.9). The deformation parameter plays the role of the coupling constant. 
Really, for small 7^0 the hopping model (1.1) expands to 



M 



H 1 = H b -± J2{bn+ibtN n + bl +1 (N n + N n+1 )b n + N n+1 b n+1 bl} + £>( 7 2 ). 



n=l 



We can consider the g-boson model as a strongly interacting boson model. On the other hand, 
the representation (1.8) suggests that the kinetic energy of (1.1) involves boson correlations, 
i.e., the hopping terms between adjacent sites depend on the occupation of those sites. These 
models are known as the strongly correlated ones. 

From this point of view the limit q — ► oo (7 — > 00) is of interest. It follows from ( [l.E| ) and 
(1.8) that the operators B, B* transform into operators 0, <ft : 



B, 



j • 



St + 1)-*, 

with the commutation relations 

[N i} <f>j] = -faSij , [Ni, (j))] = tySij , [<j>i, 4>\] = niSi. 



(1.10) 



(1.11) 



where iTj is the local vacuum projector, irj — (| 0)(0 One may verify that <fi, (fy , and can 
be expressed in terms of the Fock states, | n), as 

00 00 00 

</> = ]T \ n)(n + l\,tf = J2 \ n+l)(n\,N = Y i n\n){n \ . 

7i—0 n—0 n=0 

The introduced operator <j> is "one-sided unitary" (or isometric), i.e., 

#t = i ; 

but 



= 1-|0)(0|. 

The operators 0, $ are studied intensively in quantum optics in connection with the phase- 
operator problem (see p5|,p6[ ] and references therein). The introduction of phase variables for 
bosons was discussed in [ 27[. The relative phase of boson fields is important in the theory of 
beam splitters p8| and Joscphson junctions E9]. 



For 7 = 00, the Hamiltonian (1.4) becomes 



and 



H = ~\ Y,i4>iK+i + <f>n4>i+i - 2iV„) - fiN, 

n=l 

[H,N] = 0, 



(1.12) 



where the total number operator N is given by (1.3). This model is called the phase model. 
It should be mentioned that the case M = 2 with £l = corresponds to the phase difference 
operator considered in pa] . 
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2 The Bethe Ansatz solution of the model 

Define the L— operator for the g-bosons hopping model at lattice site n as 



e 



A 



L n (\) = ( xBn ) , (2.D 



where B n , B^ are the g-bosons ([L2]), \ = \/l — cp 2 , and A 6 C is the spectral parameter. This 
L-operator satisfies the bilinear relation: 

R(X, n)L n (\) g) L n (/i) = L n (/*) ® L„(A)fl(A, /i), (2.2) 

with the (gauge transformed) trigonometric i?-matrix 



J?(A,/i) 



/ /(m,A) \ 

<?( M ,A) g 

q- 1 ff (/i,A) 

V o o o /CM) / 



with the matrix elements being defined by the functions 



(2-3) 



/(A /i) — smn (A ~ A* + 7) . ,^ v _ sinli7 

sinh(A — p) ' sinh(A — //) 

The monodromy matrix is introduced in the usual way as 

T(A) = L M (A)...Li(A) = ( ) . (2.5) 

The 16 commutation relations of its matrix elements are then given by the i?-matrix 

R(X, /i)T(A) ® T{p) = T{fi) <g> T(X)R(X, p), (2.6) 

Let us write explicitly the relations important for deriving the Bethe Ansatz 

qA(X)B(fj,) = f(X, fx)B^)A(X) + g(ji, X)B(X)A(ji), (2.7) 

qD{X)B(p) = f(n, X)B(p)D(X) + g(X, (i)B(\)D(p), 
[B(X),B(fx)]=0, 
C(X)B(p) - q- 2 B(^C(X) = q- 1 g(X, f i){A(X)D^) - A(n)D(X)}. 
The matrix trace of the monodromy matrix is the transfer matrix 

t(A) = A(X) + D(X), (2.8) 



It follows from (2J3) that 

[r(A),r(/x)] =0. (2.9) 
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The Hamiltonian (1.1) of the g-boson hopping model is expressed by means of trace identities 
in terms of the transfer matrix 



1 — 2A 9e r(A) 



2 d\ 
and commutes with the transfer matrix 



1 2A de~ MX r(X) 

\X— — oo 2 |A— oc 



(2.10) 



[F„r(A)]=0. (2.11) 
The eigenvectors of the transfer matrix and hence of the Hamiltonian are of the form 



N 



i; N (\ 1 ,...,X N )) = l[B(\ j ) |0>, 



(2.12) 



where the vacuum state | 0) is defined in (L6), and parameters Xj satisfy the Bethe equations: 



a 2MXj _ TT f(^k,Xj) . a — i M 



(2.13) 



These states arc called Bethe vectors. One can also construct dual Bethe vectors 

N 



<oi n^cA^, 



where {Xj} satisfy the same Bethe equation (2.13). The vacuum vector |0) and the dual vacuum 
(0| are eigenvectors of the operators A(X) and D(X), with the "vacuum" eigenvalues a(A) and 
d(X), respectively. For the q-boson hopping model one has 



a(A)=e M \ d(X) = e 



-MA 



(2.14) 



The eigenvalues 9n of the transfer matrix for the eigenvectors ( 2.12 ) are 

T (^) I ^n) = N (n, {Xj}) | ip N ), 

N N 



„ N ti 



( M , {Xj}) = e M » J] f(Xj,n) + e~ M " J] A,). 

3=1 3=1 



In the explicit form the equations ( 2.13 ) are rewritten as 



JY 



JM Pj = TT Sin(£^+z 7 ) 
ll si n(^-i 7 ) ; 



where we have introduced the momenta 



p = -2iX. 



(2.15) 



(2.16) 



(2.17) 
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Later it will be convenient to use the parameters A for arbitrary sets {A} and momenta p for 
the sets satisfying the Bethe equations. In the scaling limit these equations go into the Bcthc 
equations for the Bose gas. 

Eigenenergies of the Hamiltonian (1.1), H q \ ippj) = En \ iPn), can be found from the 
equations ( |2.15| ), ([O]) and fl2.10| ) 



N 



E N = J2 h (P^ Hp) -2 sin 2 (p/2). 



(2.18) 



fe=i 



Eigenenergies of the Hamiltonian ( |1.4|) , Hh \4>n) — En \iPn), are 



N 



fc=i 



3 The phase model 



The phase model is the model defined by the Hamiltonian (1.12). It belongs to the class 
of strongly coupled correlated boson models since the hopping terms between adjacent sites 
depend on the occupation of those sites which is evident from the representation (1.1C) of the 4> 
-operators in terms of the ordinary bosons. 



The L— operator |10| of the phase model is obtained by putting j = oo in (2.1) 



L n (X) 



-ip/2 



(3.1) 



where <pln4>n are defined by ( 1.11 ). This operator satisfies the bilinear relation 

R(p, s)L n (p) (g> L n (s) = L n (s) (gi L n (p)R(p, s), 
with the 4x4 R- matrix R(p, s). The non-zero elements of the R- matrix are 

#11 (P, s) = R 4 a(p, s) = f(s,p), 
R22(p,s) = R 3 3(p,s) = g(s,p), 
R23(p,s) = 1 



and 



f(p,s) 



e 2 



2sin(^) 



g{p, s) 



2sin(^£) 



(3.2) 



This i?-matrix is obtained as the limit of the i?-matrix (2.3). 
The Bethe equations for the model are 



exp{i(M + N) Pj } = (-1) 



N-l 



exp 



{i^Pk}, 



(3.3) 



k=l 
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(j = 1, N) and are exactly solvable: 

Pj Af + A ' 1 J 

where Ij are integers or half-integers depending on N being odd or even. The Bethe vectors 
form a complete orthogonal basis. 

The A-particle eigenenergies of the Hamiltonian H — JxN ( 1.12 ) are 

N 

En = £>(Pk) ~ fi>> h(p) = 2sin 2 (p/2). (3.5) 

fe=i 

Here jj, is the chemical potential, < ft < 1. For M ~ 2 and p, = this result coincides with one 
obtained in . 

The thermodynamics of the model is handled in the standard way. It will be considered for 
the case of the zero total momentum P = ^2^ =1 Pk — 0. The state of the thermal equilibrium 
of the model at finite temperatures is determined through the solution of the nonlinear 
integral Yang- Yang equations which are drastically symplified in our case 

ln(l + e^< rt )dp, (3.6) 
2np(p)(l + e^ p >) = 1 + / p(p)dp. 



The function p(p) is a quasi-particle density while e(p) is the excitation energy. The pressure is 
then 

V = (27T/3)- 1 / ln(l + eP e &>)dp, (3.7) 



and the density is 



So we have 



dv r 

D= ^=lv {p)dp - 



(3.8) 



e{p) = h(p) -pl-V (3.9) 
and the quasi-particle density has the Fermi-like distribution 

2tt P ( P ) = (1 + D)(l + e ^)- 1 . (3.10) 

At zero temperature (/3 _1 = 0) the ground state is the Fermi zone, —A < p < A (A < n), 
filled by the particles with the negative energies eo(p). The pressure and density are now 

/A pA 

e (p)dp, D Q = / p(p)dp, (3.11) 
-A J -A 

where (see (3^) and ( 3 . 1 0| ) ) 

eoip) = h(p) -p-Vo, e (±A) - 0; (3.12) 
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2Trp (p) = (1 + Do). 
The bare Fermi momentum A is expressed as a function of density 



The Fermi velocity v is given as 



£ ° (A) -(l + Do^sin^ 



2^Po(A) v uy 1 + Da' 

If A — > [fi — > 0), then Do — > and 'Po —* as should be expected. If A — > 7r (/2 — > 1), all 
the vacancies are occupied by particles D — > oo, Vo — > 1, and the phase model becomes the 
classical AT Y chain in this limit [jl], Q . 

4 Scalar products and norms in the phase model 

Here we calculate the scalar products of the states produced by the operators B(X) and the 
norms of Bethe vectors using the standard procedure |H| |). It is necessary to calculate the 
correlation functions. Consider first the scalar products in the g-boson hopping model (|l . l[) 



§({Xf }, {Afe }) = (0|C(Af ) . . . C(\%)B(\Z) . . . B(Xf )|0>, (4.1) 

where C(A), B(\) are the matrix elements of the monodromy matrix ( |2.5[ ) and {A^} and {Aj?} 
are the sets of arbitrary spectral parameters (Bethe equations are not imposed). Using the 



commutation relations (2.7) one gets 

N N 



saxf}, {Afc }) = e n «( a «) n d (^) K N ( \ f A \ \ X X B D \ ) ■ (4.2) 

a n ™— i a— i \ ' J I ' / 



AD a=l p=l 



Here the sum is taken over all the partitions of the set {^f} U {A^} into two subsets 

{Af}u{A^} = {A^}u{A^}. 
The number of elements in each set equals N: cardjA^} = card{A^} = card{A„ } = cardjA^} 



N. Functions a(A) and c?(A) are the vacuum eigenvalues of the operators ^4(A) and -D(A) (2.14) 



{A c } {A B } 



The coefficient Kn ( j^cj {\ B } ' " n i§ nes *" coefficient) can be represented in terms 



of the partition function of the six- vertex model corresponding to the R- matrix (2.3). It was 
shown in p3fl that this partition function is just the partition function of the six- vertex model 
corresponding to the XXZ i?-matrix. The explicit expression for this partition function was 



given in 30 . This function depends on two sets of parameters {A a } and {ut}. 



N N 

f] J] sinh(A a -Vk-%) sinh(A - v k + ~j) 

?n = (-lf a = lk =i , det N M, (4.3) 

J] smh(A a - At) J] smh(ui - v k ) 

a<0 k<l 
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where the N x N matrix M has the following form 



, . sinh7 

M a ,k = — -TTT 7 s ■ wT -TT- ( 4 - 4 ) 

smh(A a - ^ fe - ^) smh(A a - f fc + £J 



The highest coefficient can be expressed as follows 



N 

\| j A cj {\ B })= e "" \ n sinh (Af - A£) ) ^({Af},{A£ + I}). (4.5) 



{A c } {A 5 } \ . 



To prove this proposition one should consider the inhomogeneous gauge transformed XXZ model 
on a lattice of N sites taking into account that the highest coefficient depends only on the R- 
matrix (see || ) . Other coefficients in (4.2) can be analogously represented as 

KN ( iS IS ) =e-^(N-n a) I ^ ^ f { ^C X c\ x 



{A } {A } 

n n fix** \ db \ \ k ( ^ ^ DC ^ \ k ( ^ ^ DB ^ 



ViGAB m£DB 



(4.6) 



Now it is not difficult to represent the scalar products for the q-boson hopping model as mean 
values of determinants depending on dual fields using the approach of |]l6|| . 

In this paper we will not, however, use this method concentrating our attention on the phase 
model only. It appears that the scalar products and correlation functions in the phase model 
can be represented without auxiliary fields. The matrix elements of the matrix Ai (4.4) the 
limit 7 — > 00 are given as 

2e 2X f 

Mj k = -^ «b +0(0, 



2AF 2A f 

e h — e 3 



JJ( e 2Af _ e 2Af) Y\ (e 2A ™ -e 2A *) 
detM = 2 N ^ ~ e2&Af +0(O, 



for the highest coefficient (4.5) we obtain 

v ( {A C } {X B } \ [ T-r, 2 \ B 2Xf \ TT / 2X C 2A°\ 



X 



k<m 



xdet N )C({\ B },{\%}), (4.7) 
exp((2iV - 1)A£ + Af ) 



'3* ~~ 2A C 2A B 

e * — e j 
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and other coefficients (4.6) are 



M fS! ) = [n«= 2X? - 2Af )n(^- 2Af )] 

i J<J k<rn 



{X A } {x } 



x (-l)(M+[Q])det„ ^({A cc }, {A AB })detjv_ no ;C({A flB }, {A AC }). (4.8) 

Here P and Q are the transpositions {X B } ->■ {A AS } U {A 155 }, {A c } -> {A- 40 } U {A DC } respec- 
tively. 

It is convenient to consider the "normalized" scalar product in the phase model 

N N 

saxf}, {A^D = <oi n c (^) n B ( A f )i°>. ( 4 - 9 ) 
k=i j=i 



where 



B(A?> = f|±, C(^, = |§. ,,10, 



From (4.2, hq) one obtains 



W?MA?}) = [ Il( e2Af " e2Af ) II -e 2A? ) ] x 

\j<l k<rn 



N 



5>l)([mKI) JJ r(A A) det„ /C({A^},{A^})det Ar _„ /C({A^},{A^}), 



x 

a.d V=i 



where r(A) = a(A)/d(A). Using the formula for the determinant of the sum of two matrices we 
can express the scalar product as a determinant of an N x N matrix 



5({Af },{A?}) = H(e^-e^) ft _ e ^) det N F({Xf }, {X c k }), (4.11) 

\j<l k<m J 

r(Af ) exp((2JV - l)Af + A?) - r{X c k ) exp((2iV - 1)A£ + Xf ) 

^i-fc = ax? ^ ■ (412) 

e j — e * 

For the phase model, r(A) = e 2XM ( |2.14j ). Hence we obtain 

e X p((2AT + 2M-l)Af +A^-exp((2AT + 2M-l)A^ + Af) 

= 2X? ^ ' (4 ' 13) 

e j — e * 

There is an important case when {X B } = {A c }. The diagonal elements of the matrix F are then 

Fjj = (M + 7V_i) e 2(M+w-i)A j5 
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Turn now to the norms of Bethe eigenvectors 

N N 

\iP( Pl ,...,p N )) N = l[B(\ ] )\0), (V(pi,...,PJv)| = (0|n c ( A i). Pi = -2iAj, (4.14) 



where {p} is a solution of the Bethe equations (3. 3). Now we have 

JV 

^• = (M + iV-l)e^-^, P = J>, 

i=i 

It is not difficult to calculate the determinant 

det^ = e^-^MCM + N) N ~ 1 , 
and for the norm of any Bethe eigenstate we obtain the formula of the Gaudin type 



7V 2 (pi, ...,Pn) = (ip(pi, ■ ■ ■ ,Pn)U>(Pi, ■ ■ ■ ,Pn)) = f JJ eip . \ eipk j 



e^'^MiM + Nf- 1 . 

(4.15) 



5 The darkness formation probability 

Here we will calculate the darkness formation probability, i.e., the probability of the states with 
no particles on the first m sites of the lattice. Formally it can be defined as the normalized mean 
value 

r(m, {p}) = Af~ 2 (pi, ■ ■ ■ ,PN){i>(pt, ■ ■ ■ ,PN)\exp{aQ(m)}\ip{pi, . . . ,PN))\ a= _ QO , (5.1) 



where Q(m) = SJLi Nj is the number of particles operator on the first m sites. 
The monodromy matrix can be represented in the following form 



T(A)=T 2 (A)T 1 (A), 

T 2 (A) =L M (X)...L m+1 (X) =(q 

Ai(A) Bi(\) 



A 2 (\) B 2 (A) 
' 2 (A) D 2 (X) 

T 1 (X)=L m (X)...L 1 (X)= ( ^ (A) Di(A) 
The bare vacuum (|l.6| ) can be represented as 

|0) = |0) n ®|0)i, 



(5.2) 
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where |0)i and |0)n are the Fock vacua for the monodromy matrices T\ (A) andT2(A), respectively. 
We denote r\(X) and r 2 (X) the ratios of the corresponding vacuum eigenvalues 

r x (A) = a 1 (X)/d 1 (X) = cxp{2mA} ,r 2 (A) = a 2 (X) / d 2 (X) = exp{2(A/ - m)A}. 

The numerator of (|5 . lh is a particular value of the following matrix element 



N N 

T({Xf}, {A£}) = <0| J] C(A£) exp{aQ(m)} J[ B(Xf )|0)| a= _ O o ) (5.3) 
fe=i i=i 



where the operators B(A) and C(A) are defined by ( 4.1C ) and {A^}, {Aj?} are arbitrary sets of 
N spectral parameters. Using the decomposition (5.2) one obtains 

N N N 

T({\f},{\%}) = l[r 1 (Xf)u(0\ II C a (Af) JjB a (Af )|0>n. 

3=1 k=l 3 = 1 



The scalar product here can be calculated by means of ( 4.11 ) , with the result 

T({Xf}AXZ})=(u(r 2Xf -e 2Xr )U(e 2XCm -^)) detr({Af},{A?}), (5.4) 

\j<l k<m J 

the matrix elements of the N x N matrix T being equal to 

x (exp((27V + 2M - 1)A£ + Af ) - exp((2iV + 2M - 2m - l)Xf + (2m + 1)A£)) . 

If the sets coincide, {X B } = {X c } = {A}, then the diagonal elements should be understood in 
the sense of the l'Hopital rule 

T n = (M + N-7U- l) e 2(M+iV-l)A j _ 

If the spectral parameters satisfy the Bethc equations the matrix elements Tjk can be represented 
in the form (pj = —2iXj and P is the sum of momenta) 

Tjj = (M + N-m- l)e i( - p - p *\ 
sin 5 (pj -p k ) 

Thus one obtains the following result for the darkness formation probability 

r(m,{p}) = T ^( } {p}) }) = (1 + D)detf{m > {p}) > (5 - 5) 
- 11 sm^±l( Pj -p fc ) 

lift = Oj.fc-— . \ — , (5.6) 

M l + D smifpj -p k ) 

where D = JV/M. 
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6 Time-dependent correlation functions 

In this section a determinant representation for the two-point time-dependent correlation func- 
tion of fields on the finite lattice is obtained. The derivation is somewhat similar to the derivation 
in the case of the XXO chain jljj] . To this end one needs first a representation for the form factors 

GN(m,{p},{q}) = — -r —j- — , (6.1) 

N(qi, . . ,,q N -i)N{pi, . . . ,p N ) 

<#>(■»■ w. M> - • -""g" ■ ■ -/f ■ M 

■A/ (gi, . . . , gw+ijA/ (pi, . . . ,pn) 
where |^({p})) are Bethe vectors and Af({p}) are their norms (4.15) 

AA(W) = W{p})M{p})> 1/a . 

The determinant representations for the form factors of the phase model were obtained in 
|p3[ (a brief derivation of these results is given in Appendix A) : 

——l 

G N (m,{ P hU}) = M- 1 (M + N)- N +iZe- ime ( -^±^ j 2 x 

x (l + ^)det JV -i( J Di-a; J D 2 )| x= o ! (6.3) 

where is a complex number depending on P and Q, \Z\ = IThe explicit form of Z is not 
written here since since it appears to be of no importance for calculating the correltors. We use 
convenient notations: 

e = Q-P, Q = X>, P = Y,Pr 

j j 

Matrix elements of matrices D\ and D% are 

D hk = e^ C S \f- Pi) + ^\{Q-p,)), (6.4) 
tan ^iPj - Qk) 2 

D 2jk = e*»i ™\f- p »\ + *nhe-p N )). (6.5) 



Analogously, for the form factor (S.2) we have (the star denotes the complex conjugation) 

x (1 + — )detjv(£>i +) - xD ( 2 +) )\ x=0 , (6.6) 

<» = e-i» + .(^l^±H+sini(e- w+I )). (6.8) 
3 tan 2(^+1 -Pj) 2 



where 
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Consider now the normalized mean value of the time-dependent product of two phase oper- 
ators 

JV JV 

<oi n c( Pfc )0 m+1 (t)4(o) n B(p fc )|o> 

/+(M,m,i) = , (6.9) 



JV JV 

n c(p fc ) n B(p fc )io) 
fc=i fc=i 



where 



,(*) = exp[ii?i]0 m exp[-zfft], 



with the Hamiltonian ii (1.12). Using the formulae ( |6 .3| ) , ( |6.6| ) we can represent this correlator 
as follows 



/+(M, m, t) = M-\M + N + 1)-^ ( M ^^ 



JV-1 



JV+1 



/V 



JV+1 



k=l 



3 = 1 



3 = 1 



k=l 



(6.10) 



(1 + — )det N (D[ + \{p}, {q}) xDi +) ({p}, {q}))\ x=0 



where e(p) is the energ y of quasiparticle e(p) = 2 sin 2 | — /i (3.5). 

The summation in (|6. 10 ) is taken over all the solutions {q} of the Bethe equations ( |3.3| ) such 
that cardjq} = N + 1. It can be seen from ( |3.4|) that Q = ^2 qj — where if is an integer, 



M 



< K < 



A I 



Then the sum over all the solutions of the Bethe equation can be rewritten 



as the sum over all such Q and the sum over all the sets of N + 1 different qk satisfying the 
following conditions 



9fe 



27T/fc + Q 

M + N + l' 



N+l 

7T < Qk < 7T, Q = ^2 QK- 

k=l 



Taking into account that the form factor is an antisymmetric function of momenta {q} one can 
make the following substitution 



(M + N 



M+N / JV+1 \ 

A ; Q 1=0 qi gjv+i \ fc=l / 



Q 1=0 qi 9JV + 1 

The determinants in ( |6.10|) can be then rewritten as follows 



(1 + ±)det N (D{ + \{ P }, {q}) - xDi + \{p}, {q}))\ x=0 



N 



E( _ 1)[s] n e ^ a rf-^\ + sin i (0 _ gQj)x 

o „-i tan 5 (<7g a -p Q ) 2 
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x(l + ±.)det N (D[ + \{p}, {q}) - xDi + \{p}, {g}))| x=0 , 

where the sum is taken over all the permutations Q of {qi, . . . , qjv+i}- We can perform this 
summation using again the antisymmetry of the form factor 

{N+l)l*g ^ tan£(g Qa -jJ ) 2 

n^ a , COS 5( -9a) . 1- ^ 

( tanl( ga - Pa ) +Sln 2( e -^- 

Taking into account that det(Di — xD?) is a linear function of x it is possible to rewrite 
( fUOl ) as 



/jv(W. m, t) = f M +f + M " Y M + % + 1 £ e^(W. *. 6), (6.11) 



AT-l M+N 



V M + N J ^ M 2 



frwCM.M,©) = exp ( -U^Pj +it^2e(pj) 

3=1 3=1 



H • ' ■ H ( fl .f ffltp[^«Ar+i - **e(«jyr+i)] + S- ) dfitW(aj)| B=0 , (6.12) 



x 

qi gw+i 



M + iV + 1 " L " T V "'" ri/J &c 
where the N x N matrix W^x) is defined as 



W(z)=W (1 ' - W^ 2) , (6.13) 



sinke-^)], (6.14) 



x / cos |(e - gg) ^ ^ i ^ 

\tan±(<j a -p 6 ) 2 

^ } = (M + TV + l) 2 exp ^ + \) qa ~ its (^ exp[i(Z - ^)q N +i - ite(q N+1 )] x 

cos±(9-g a ) . 1 \ / cos±(9 - gjy+i) .1 A 

f- + sin -(0 - g„) J I f- + sin -(9 - q N+1 ) . (6.15) 

tan 2(qa-Pa) 1 J \tan 5 (qjv+i - Pb) 2 / 

The matrix is of rank one and hence detjvW(x) is a linear function of x. 
Let us introduce the following functions 

= M ~ 1 X! ex P^g ~ »fe(g)]> (6-16) 
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l(l,t,p) = 



M 



i V 

■N + 1 ^ 



exp[iZg — its(q)) 
tan i(g — p) 



(6.17) 



d{l,t,p) = 



c^q>[ilq — ite(q)] 
(M + N + 1) 2 ^ sin 2 i((z-p) ' 



E 



(6.18) 



The summation in these formulae is made over all the permitted values of q for fixed, thus 
these three functions depend also on 6. For t = they can be calculated explicitly. 



9(1,0) = Sifi, 



e(l, 0,p) = i(l - 8 L0 - i tan -(p - 6))e^, 
d(l,0,p) = cos- 2 i(0 -p)e«* + A/+ ^ +1 ^(Z,0, P ). 



(6.19) 
(6.20) 
(6.21) 



Using these functions we can perform the summations in (6.12), bearing in mind also the fol- 
lowing relation 



cot \(l - Pa) cot hq- Pb )= cot i(p - p b ) 



1/ \ 1, 

cot o W ~ Pa) ~ C0t o ' 9 ~~ W 



- 1. 



Performing the summation in each row and taking into account that detW(x) is a linear function 
of x we have 

h+({p}, I, t, 9) = (g(m, t) + dct N [S - xR+] | a=0 . 



(6.22) 



The matrices S 1 and i? + are given as 



S, 



1 



ah 



i 



M + iV+l itani(p a -p b ) ^2 



:(et(l,t,p a ,Q) + e + (l,t,p a , Q))e-(l,t,p b )- 



- i (e+(l,t,p b ,e)+e + (l,t,p b ,Q))e_(l,t,p a ) 



-g(l,t)e-(l,t,p a )e-(l,t,p b ) + -(e£(l,t,p a ,Q) - e + (l,t,p a ,Q))e^(l,t,p b )+ 
+ ^(e+G:*, ?>&,©) - e+(l,t,p b ,e))e-(l,t,p a ) } , for a ^ 6, 



(6.23) 



1 



1 



5 aa = - ( d(l,t, Pa ,e) + ^(e~ ze d(l + l,t, Pa ,e) + e™d(l - l,t, Pa ,Q)) I x 



xe-(l,t,p a )e-(l,t,p a ) + 



1 



M + N+l 



{-g(l,t)e-(l,t,p a )e^(l,t,p a )+ 



+ i(ef(l,t,p a ,Q) - e + (l,t,p a) Q))e_(l,t,p a )}, 



Kb = M + N + 1 e +( l i t 'Pa,Q)e + {l,t,p b ,Q). 



(6.24) 
(6.25) 
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The functions e_(Z, t,p), ei(i, t,p, Q), e , (I, t,p, 0) are defined as 



e-(M,p) = exp I -4p + 4 £ (p) ) ' 



(l,t,p,e) = -e^(l,t,p)x 



x ((e(/, t,p) + e~ ie e(l + 1, - t) - e' 40 ^ + 1, f))) , 
e+(l,t,p,9) = -e-(l,t,p)x 



(6.26) 



(6.27) 



x ((e(Z, t,p) + e l6 e(Z - 1, t,p)) + i(g(l, t) - e lB g{l - l,t))) 
In the equal-time case one has explicit expressions for these functions 



(6.28) 

(6.29) 

(6.30) 
(6.31) 

Equations (|s7Tl| ) and fl6.22H6.25l) give a determinant representation for the correlation func- 
tion ( |6.9| ) on the finite lattice. 

The calculation of the two-point time-dependent correlation function 



e_(I,0,p) = e- 1 **' l 
e+(Z,0,P,e) = z(l-Me 1 ^, 



N 



(oi n c(pk 



/jv(M. m '*) 



/c=i 



^t)M0) ft B(p fc )|0> 



TV iV 

(oi n c(p fc ) n B(p fc )io) 

fc=l fc=l 



(6.32) 



is quite similar and we give only the final result 

JV-1 M+N-l 



/jv(to' TO ! f ) 



M + 7V- 1 
M + iV 



£ - E eiM9 ^(W.M,6), (6.33) 

i=i e 



fc^({p}, I, t, 6) = ^det w [5 + ariZ - ] 



x=0j 



D— 

afc ~ M + N - 1 



e-(l,t,p a )e^(l,t,p b ). 



(6.34) 
(6.35) 



One shoul d note that in this case the fac tor M + N + 1 entering the definitions of the matrix S 
( |6.23 6.24 ) and the functions e, d and g ( |6. 16 - 6. 18 ) should be replaced by M + N — 1. 
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7 The thermodynamic limit 



Let us consider the thermodynamic limit (the total number of sites M — * oo, number of particles 
N — > oo, the density D = N/M remains finite) of the correlation functions obtained in the 
previous sections. In this limit one should replace all the sums over momenta by the integrals 
and the determinants oi N x N matrices by the Fredholm determinants of the corresponding 
integral operators acting on the functions on the inter val [— tt, it] H. 



The result for the darkness formation probability (5.5) is 

r(m) = (l + D)det(I-T), (7.1) 
where I is the identity operator and T is an integral operator 

Cf/)(p)= / T(p,q)f(q)dq, 



with the kernel 

1 sin m±l(p-q) 

Tip, q) = . r-p(q)- (7.2) 

1 + D sm i (p — q) 



Here p(q) is the quasi-particle Fermi-like distribution function (3.1C) 



Pip) - ^(1 +D)(1 + expipeip)))- 1 , (7.3) 



where is the temperature, the total density D is defined by (3J3) and the energy e(p) 
is the solution of the non- linear integral equation (|3.6| ). After symmetrizing the kernel the 
representation for the darkness formation probability can be rewritten as follows 

r(m,/3) = (1+D)det(/-M), (7.4) 

where M is an integral operator with the kernel 

m+l , 



1 rr St sin2^i(p-g) 



? ) = — v>P)— VVP), (7-5) 
Zn sm — gj 

and v(p, 0) = (1 + exp(/3e(p))) _1 is the Fermi weight. 

At zero temperature (/3 _1 = 0) the Fermi weight becomes the step function equal to zero 
outside the Fermi zone and equal to one inside it. Thus one has 

r (m) = (1 + D )det(7 - M ), (7.6) 

where Mo is an integral operator 

A 

(M /)(p)= / M {p,q)f{q)dq, (7.7) 

-A 
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with the kernel 



M (p, q) 



1 sm^±l(p-q) 
2-7T sini(p — q) 



(7.8) 



and the Fermi momentum A is defined by the equations ( 3.12| ). 

Consider now the two-point time-dependent correlation function (6.9) in the thermodynamic 
limit. The equation (3.11) takes the form 



/( ± )Kt,/3) = exp(±-^-)l±^f; / e ^-^h^(i,t,p,e)de 

The functions {I, t, 0, 9) can be written as Fredholm determinants 

ftW {I, t, p, 9) = (g(1, t) + ^\ det(J + V - xR + )\ x = , 

h^(l,t,0,@) = ^-dct(i + V + xR-)\ x=0 , 



where V and R are integral operators 



(7.9) 



(7.10) 
(7.11) 



(Vf)(p) 



1 

27 



V(p,q)f(q)dq, 



(R ± f)(p) = ^ I R ± (p,q)f(q)dq, 



with kernels 

V(p,q) 



tan 5 (p — g) V 2 



1 



- t,p, /3, 6) + E+ (I t, p, (3, G))E_ (I, t, q, /?)- 



(7.12) 



- - (E+ (I, t, q, (3, 0) +E+(l, t, q, (3, 9))£_ {I, t, p, f3) ) - 
-G(l,t)E_(l,t,p,(3)E_(l,t,q,0) + 7 -(E+(l,t,p,(3,0)-E+(l,t, P ,P i e))E_(l,t,q,p) + 



+ -(E+(l,t,q,(3,e)-E-(l,t,q,(3,e))E_(l,t,p,(3), 



R+(p, q) = E+(l, t,p, 0, Q)E7(l, t, q, (3, 9), 



(7.13) 
(7.14) 

R~(p, q) = E_ (I, t,p, 0)E- (I, t, q, (3). (7.15) 
The functions G(l,t), E-(l,t,p, (3), E~\_{l,t,p, (3,0) and ET.(l,t,p, (3,@) are defined as follows 



G(l,t) 



2tt 



exp(ilq — ite(q))dq, 



(7.16) 
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7T 

E(l,t,p,Q) = ±w.p. f ex P(^g ~ »te(g)) dq + tan l(p _ 0) exp(aZp - ife(p)), (7.17) 
27r 7 tan^g-p) 2 

— 7T 

£_(Z, /?) = ^(p,/9)ex P ((-i^p + i^e(p)) ), (7.18) 

£?+(/, t,p, J 9,e) = ^_(/,t,p, ) 9)x 
x ((E(l, t,p, 9) + e- 4e S(? + 1, t,p, 9)) - i(G(l, t) - e- ie G(l + 1, t))) , (7.19) 
E+(i,t,P,/?,e) = ^_(Z,t,p,/3)x 

x ((£(Z,t,p,9) + e* e £(Z - l,t,p,e)) + i(G(/,t) - e m G(l - l,t))) . (7.20) 

One should note that the function E(l,t,p,&) is singular at the points p = 6 ± ir but the 
functions E+(l,t,p, (3,Q) and E+(l,t,p, /?, 9) entering the kernels have no singularities being 
well defined for all 6 and p. 

In the case of equal-time correlators these functions can be calculated explicitly 

G(i,0) = <y,, (7.21) 

£L(Z,0,p,/3) = V^P^e"'*", (7.22) 



E+(l,0,p,p,S) = i(l - 5,, )\Ap)e i ^, (7.23) 

(Z, 0,p, /?, 9) = i(l - t^)^ 6 ^ y/Ufafle'b, (7.24) 

and the kernels V(p, q) and i?(p, g) are polynomials in e 10 . Hence, h^\l, 0, (3, 9) can be repre- 
sented as Taylor series in e ie and 

s i(ro - i)e /i (±) (Z,0,/3,e)(ie = for m>Z, (7.25) 



so that 

f( ± \m,O,0)=e X p(± T ^\ /e i(m -° e /i (±) («,0,/?,e)de. (7.26) 

The functions ftW (Z, 0, /?, 9) can be represented as determinants of simple operators: 

h^\l, 0, /3, 9) = ^det(J - fi + zf*)^, (7.27) 
where the kernels of the operators t) and f 1 * 1 are 

Gin I — — 



, , ^ sin i±l(p - g) + exp[z(9 - £±2)] sin ^(p - g) _ 

v(p, g) = vMp) —r ( r VK?)> (7-28) 

sm 5 (p- (?) 
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r + (v,q) = V^y (e ' 9) ^ {p+q) yMq), (7.29) 

r- [p, q) = VHP)e~ lHp+q) yfifij. (7-30) 

Again at zero temperature the Fermi weight v{p, 0) becomes the step function and the correlators 
have the following form 

/r ) (m,0) = exp(± T ^-ji±^g / e^ e ^(i > 0,e)d8 ) (7.31) 



h( +) (l,t,Q)= (G(l, t) + ^ det(I + % - xR+)\ x = , (7.32) 
h { -\l,t,e) = — det{i + V +xRo)\ x=0 , (7.33) 



where Vo and R^ are integral operators 



A 

1 



{V f)(p) = — / V (p,q)f(q)dq, 



-A 
A 



(R±f)(p) = ±- J R±(p,q)f(q)dq. (7.34) 

-A 

The kernels of these integral operators acting on the interval [—A, A] are given by the equations 



( 7.13 - 7.15 ) after putting formally v(p, /3) = 1. 



Conclusion 

In this paper we have represented the correlation functions for the phase model as determinants. 
In the thermodynamic limit these are the Fredholm determinants of " integrable integral oper- 
ators" pS3fl . These representations will allow us to derive classical integrable equations for the 
correlators and to evaluate their large time and distance asymptotics. 

The model considered is not a free fermion model but the correlation functions are represented 
in an explicit form not involving auxiliary dual fields. It is interesting to mention that the 
corresponding point exists also for the XXZ Heisenberg chain with an infinite anisotropy. It 
is natural that our technique can be applied also at this point, and the expressions for the 
correlators will be considerably simpler than in the general case. We are going to present this 
results somewhere. 
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8 Appendix: Form factors 

In this Appendix we calculate the matrix element 



N-l N 

G({Af},{Af}) = (0| J] C(A^)^nB(Af)|0). 

3=1 



(A.l) 



k=l 



Using (5.2) with m = M — 1 and taking into account that in this case B-i = (/>Jj, C2 — <Pm we 
obtain the representation 



G({Af},{A^}) = ^n^(Af)i(0|n c i( A f)Il B i( A f)l ) 



lX 



I, II I 



n^nxyWnxAcrinn 



1 



11 d 2 (Ag) jj d 2 (Ag) y j jj e 2 *? - e 2 ^ 



nn 

V I II 



„2A, C 



(A.2) 



where the sum is taken over all the partitions of the sets {A B } and {A } into two subsets {Af}, 
{A^} and {Ap}, {Ag}, respectively, satisfying the following conditions 



card{Af } = card{Af } = n±, 



m = o,i,...,JV-i, 



card{An} = card{Ag} + 1 = N - ni = n 2 . 



Using (4.11) we can rewrite (A.2) 
G({Af},{AP})= 



e 2X r - e 2y ' 



')TJ( e ^- e ^) 



\3<l 



k<m 



N N-l 

exp(^Af + £ A£)£(-l)«™nexp(2(M + rH - l)Ag)x 
j=i fc=i 1, 11 11 



n 

3<i 



' 2A n . 2A^ > 

e J — e u i 



JJ (e 2AS, _ e 2An m ) det „ 1 ^ 1 ({Af},{Af}), 



ra<k 



where 



exp(2(Af + M - l)Af ) - cxp(2(^ - m + l)Af + 2(M + m - 2)Ap) 



2A fl 2A C 

e ^ — e * 



Let us consider the function in the left hand side of (|A.3|) 



(A.3) 



(A.4) 



2A T C - 



2A 



G 2 ({A B },{Ag})=n ex p( 2 ( M +^- i ) A ii)ii( e j ■ -« 2A "0 n( e "^ - e 

II j</ rn<k 

These products can be rewritten as determinants 

Ga({Ag},{Ag}) = JJejcp(2(M + m - l)Ag)det w _ ni M{Ag})det w _ ni _i$({Ag}), 
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where the matrices g and h are given by the formulae 

It is convenient to introduce the [N — ni) x [N — ni) matrix h 

hik = h kl = S lk , hjk = g 3 -i,k-i for j > 1, k > I, 
det N - ni h = detTv-m-15- 

Then it is easy to see that 

G 2 ({Ag}, {Ag}) = det^^H^X*}, {Ag}), 
H 2jl = 1, 

rap(2(m +M- l)A^_j + 2(N - n x )Xf ) - cxp(2Af + 2(M + N- 2)\°_ 1 ) 
^2 Jfc = ^ ^ , (A.5) 

e j — e fc - 1 



for k > 1. The representation ( A. 3 ) can be rewritten now in the following form 

G({Xf } ,{A^}) = ( U(e^f ) [] (e 2 ^ _ * 

\j<l fc<m 

iV iV-1 

exp(£ Af + £ X^i-l^^detMiXfjAX?}^ 
3 '=i /c=i i,n 

xdet W -mff2({Ag},{Ag}). (A.6) 

To use the Laplace formula for the determinant of the sum of two matrices it is necessary to 
introduce four dual fields [|l6): ipf, ipj, j = 1, 2, 3, 4 

[V'j » ^k] = ' ^ ] = o, , V'fc ] = <5j,fc . 

acting in the dual Fock space with the vacuum |0) 

^|0)=0, (0|^+ = 0. 
Using the commutation relation for the dual fields it is easy to prove the following relations 
det ni H 1 ({Xf},{Xf}) = (0|det ni £i({Af},{A?})|0), 
deW_ ni F 2 ({Af I }, {Ag}) = (0|detjv_ ni e 2 ({A£}, {Ag})|0), (A.7) 

-(exp(2(JV + M- l)Af)- 



_ CXp(-0^ +-0l) 



\B i If/If o^C | C\C 



exp(2(iV + l)Af + 2(M - 2)A^T + 2(A^ - Xf + ^)), (A., 
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= G " P i^ 4+ tt 3) ex P (2(7V + M - 1)A£_ X + 2(Af - A£L X )(^ + <M)- 



2A B 2Af , 

e j — e 



exp(2Af + 2(M + iV-2)A^_ 1 ), for fc > 1. 



(A.9) 



Now ( A.6 ) can be rewritten in the form 
G({Xf},{^})= |II( e 



2X" 2\: 
j _ e i 



A? 



N-l 



k<rn 



2\t „2A 



e^) 



exp(^Af + ^ \C)(0\det N g({\ B },{\ c })\0), 

3 = 1 



(A.10) 



fe=i 



9jl = 1, Gjk = Qi 3 , h ^ + Qi ik for k > 1. 



The determinant in (A.1C) is represented as a sum of minors 



N 



N-l 



(0|det^({A B },{A c })|0) =5>xp(£Af - ^ A fc) x 

;=i jj^i fe=i 

x (-l)'(0|det JV _ 1 F({Af,...Af_ 1 ,Af +1 ,...,Af}{A c })|0). (A.ll) 
The mean value of the determinant of the matrix V can be expressed without dual fields, 



(0|det^7({Af},{A5r})|0) = det^({Af},{ASr}). 
Now we can use this relation to rewrite G({X B }, {A^ }) also without dual fields: 
(0\det N g({X B }, {X c })\0) = det N H({\f}, {A^}), 



(A.12) 



H n = 1, 



exp(2(7V + M - l)Af ) - exp(2Af + 2(M + N- 2)\%_ 1 ) 



2A 2Af , 

e 3 — e fc - 1 



for k > 1. 



(A.13) 



Finally we have 



G({Af }, {A^}) = ]> 2Af - ^ ) II - e ^ 



\3<l 



k<ra 



N 



N-l 



expQT Af + J] X%)det N H({Xf}, {A^}). 



fe=i 



(A.14) 
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If {A^}, {Aj?} are solutions of the Bethe equations (|3.3[) the matrix H can be rewritten as 



p i{P-Pj) 4- pi{Q-qk-i+Pj) 
H » = e^-e^ ' (A ' 15) 

where 

TV AT-1 

Pj = -2iXf, q k = -2i\%, P = J2pj, Q=Y,1k- 

3 = 1 k=l 



Using (4.15) we obtain the normalized form factor 



G N (M, M, {q}) = Z M M-\M + N)- N +i Lji±^-S\ 2 detH({p}, {q}), (A.16) 

where Zm is a complex number depending on P 1 Q and M, \Zm\ = 1- 

If {A^},{A^} are solutions of the Bethe equation then, using the shift operator, it is also 
possible to calculate the normalized matrix elements of any operator </> m . 

G N (m,{p},{q}) = Z m M-\M + N)- N+ i J ' detH({p},{q}), (A.17) 



Z rn = e <M-m){Q-P) z 



M ■ 



Notice that this representation is equivalent to the formulae ( p . 3[ ) , ( J6 . 4[ ) . It is evident that now 
one can also obtain a representation for the normalized matrix elements of the operators ^J^, 
using the following relation 

G { x>(m,{p},{q}) = G* N+1 (m,{q},{p}), 
where star means the complex conjugation. Thus one gets for this form factors the representation 
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